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Unusually High Thermal Conductivity of Carbon Nanotubes 
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Combining equilibrium and non-equilibrium molecular dynamics simulations with accurate carbon 
potentials, we determine the thermal conductivity A of carbon nanotubes and its dependence on 
temperature. Our results suggest an unusually high value A«6, 600 W/m-K for an isolated (10, 10) 
nanotube at room temperature, comparable to the thermal conductivity of a hypothetical isolated 
graphene monolayer or diamond. Our results suggest that these high values of A are associated with 
the large phonon mean free paths in these systems; substantially lower values are predicted and 
observed for the basal plane of bulk graphite. 
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With the continually decreasing size of electronic and 
micromechanical devices, there is an increasing interest 
in materials that conduct heat efficiently, thus prevent- 
ing structural damage. The stiff sp 3 bonds, resulting in a 
high speed of sound, make monocrystalline diamond one 
of the best thermal conductors M. An unusually high 
thermal conductance should also be expected in carbon 
nanotubes , which are held together by even stronger 
sp 2 bonds. These systems, consisting of seamless and 
atomically perfect graphitic cylinders few nanometers in 
diameter, are self-supporting. The rigidity of these sys- 
tems, combined with virtual absence of atomic defects or 
coupling to soft phonon modes of the embedding medium, 
should make isolated nanotubes very good candidates for 
efficient thermal conductors. This conjecture has been 
confirmed by experimental data that are consistent with 
a very high thermal conductivity for nanotubes [Q. 

In the following, we will present results of molecular 
dynamics simulations using the Tersoff potential || , aug- 
mented by Van der Waals interactions in graphite, for the 
temperature dependence of the thermal conductivity of 
nanotubes and other carbon allotropes. We will show 
that isolated nanotubes are at least as good heat con- 
ductors as high-purity diamond. Our comparison with 
graphitic carbon shows that inter-layer coupling reduces 
thermal conductivity of graphite within the basal plane 
by one order of magnitude with respect to the nanotube 
value which lies close to that for a hypothetical isolated 
graphene monolayer. 

The thermal conductivity A of a solid along a partic- 
ular direction, taken here as the z axis, is related to the 
heat flowing down a long rod with a temperature gradient 
dT/dz by 
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where dQ is the energy transmitted across the area A in 
the time interval dt. In solids where the phonon contri- 
bution to the heat conductance dominates, A is propor- 
tional to Cvl, the product of the heat capacity per unit 



volume C, the speed of sound v, and the phonon mean 
free path I. The latter quantity is limited by scattering 
from sample boundaries (related to grain sizes), point 
defects, and by umklapp processes. In the experiment, 
the strong dependence of the thermal conductivity A on 
I translates into an unusual sensitivity to isotopic and 
other atomic defects. This is best illustrated by the re- 
ported thermal conductivity values in the basal plane of 
graphite || which scatter by nearly two orders of mag- 
nitude. As similar uncertainties may be associated with 
thermal conductivity measurements in "mats" of nan- 
otubes we decided to determine this quantity using 
molecular dynamics simulations. 

The first approach used to calculate A was based on 
a direct molecular dynamics simulation. Heat exchange 
with a periodic array of hot and cold regions along the 
nanotube has been achieved by velocity rescaling, follow- 
ing a method that had been successfully applied to the 
thermal conductivity of glasses 0. Unlike glasses, how- 
ever, nanotubes exhibit an unusually high degree of long- 
range order over hundreds of nanometers. The perturba- 
tions imposed by the heat transfer reduce the effective 
phonon mean free path to below the unit cell size. We 
found it hard to achieve convergence, since the phonon 
mean free path in nanotubes is significantly larger than 
unit cell sizes tractable in molecular dynamics simula- 
tions. 

As an alternate approach to determine the thermal 
conductivity, we used equilibrium molecular dynamics 
simulations ||[| based on the Green-Kubo expression 
that relates this quantity to the integral over time t of 
the heat flux autocorrelation function by |l(| 
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Here, ks is the Boltzmann constant, V is the volume, T 
the temperature of the sample, and the angled brackets 
denote an ensemble average. The heat flux vector J(t) is 
defined by 
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dt 

where Ae^ = e^— < e > is the excess energy of atom i 
with respect to the average energy per atom < e >. i"; 
is the position and Vj the velocity of atom i, and r.y = 
rj — rj. Assuming that the total potential energy U — 
J2i u i can be expressed as a sum of binding energies Uj 
of individual atoms, then fy = — VjUj, where Vj is the 
gradient with respect to the position of atom i. 

In low-dimensional systems, such as nanotubes or 
graphene monolayers, we infer the volume from the way 
how these systems pack in space (nanotubes form bundles 
and graphite a layered structure, both with an inter-wall 
separation of «3.4 A) in order to convert thermal conduc- 
tance of a system to thermal conductivity of a material. 

Once J(t) is known, the thermal conductivity can be 
calculated using Eq. (|^). We found, however, that these 
results depend sensitively on the initial conditions of each 
simulation, thus necessitating a large ensemble of simu- 
lations. This high computational demand was further 
increased by the slow convergence of the autocorrelation 
function, requiring long integration time periods. 

These disadvantages have been shown to be strongly 
reduced in an alternate approach 11 that combines the 
Green-Kubo formula with nonequilibrium thermodynam- 
ics Jl^,[l3| in a computationally efficient manner In 
this approach, the thermal conductivity along the z axis 
is given by 
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where T is the temperature of the sample, regulated by 
a Nose-Hoover thermostat Jig ]) and V is the volume of 
the sample. J z (F e ,t) is the z component of the heat 
flux vector for a particular time t. F e is a small fictitious 
"thermal force" (with a dimension of inverse length) that 
is applied to individual atoms. This fictitious force F e 
and the Nose-Hoover thermostat impose an additional 
force AFi on each atom i. This additional force modifies 
the gradient of the potential energy and is given by 



AFj = Ae^Fg — fyTryFe) 
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Here, a is the Nose-Hoover thermostat multiplier acting 
on the momentum pi of atom i. a is calculated using the 
time integral of the difference between the instantaneous 
kinetic temperature T of the system and the heat bath 
temperature T eg , from a = (T — T eq )/Q, where Q is the 
thermal inertia. The third term in Eq. (j^) guarantees 
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FIG. 1. (a) Time dependence of the axial heat flux J z [f) in 
a (10, 10) carbon nanotube. Results of nonequilibrium molec- 
ular dynamics simulation at a fixed applied thermal force 
F e = 0.2 A -1 , are shown at temperatures T = 50 K (dashed 
line), 100 K (solid line), and 300 K (dotted line), (b) Time 
dependence of J z (t)/T, a key quantity for the calculation of 
the thermal conductivity, for F e = 0.2 A" 1 and the same tem- 
perature values, (c) Dependence of the heat transport on the 
applied heat force F e in the simulations for T = 100 K. The 
dashed line represents an analytical expression that is used 
to determine the thermal conductivity A by extrapolating the 
simulation data points A for F e — >0. 
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FIG. 2. Temperature dependence of the thermal conduc- 
tivity A for a (10, 10) carbon nanotube for temperatures below 
400 K. 



that the net force acting on the entire iV-atom system 
vanishes. 

In Fig. |l| we present the results of our nonequilibrium 
molecular dynamics simulations for the thermal conduc- 
tance of an isolated (10, 10) nanotube aligned along the z 
axis. In our calculation, we consider 400 atoms per unit 
cell, and use periodic boundary conditions. Each molecu- 
lar dynamics simulation run consists of 50,000 time steps 
of 5.0xl0~ 16 s. Our results for the time dependence of 
the heat current for the particular value F e — 0.2 A -1 , 
shown in Fig. 0(a), suggest that J z (t) converges within 
the first few picoseconds to its limiting value for t^oo in 
the temperatures range below 400 K. The same is true for 
the quantity J z (t)/T, shown in Fig. |l](b), the average of 
which is proportional to the thermal conductivity A ac- 
cording to Eq. (|J). Our molecular dynamics simulations 
have been performed for a total time length of 25 ps to 
represent well the long-time behavior. 

In Fig. |l](c) we show the dependence of the quantity 
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on F e . We have found that direct calculations of A for 
very small thermal forces carry a substantial error, as 
they require a division of two very small numbers in 
Eq. (JsJ) . We base our calculations of the thermal conduc- 
tivity at each temperature on 16 simulation runs, with F e 
values ranging from 0.4—0.05 A -1 . As shown in Fig. |l|(c), 
data for A can be extrapolated analytically for F e ^0 to 
yield the thermal conductivity A, shown in Fig. ^|. 

Our results for the temperature dependence of the 
thermal conductivity of an isolated (10, 10) carbon nan- 
otube, shown in Fig. ||, reflect the fact that A is propor- 
tional to the heat capacity C and the phonon mean free 
path At low temperatures, I is nearly constant, and 
the temperature dependence of A follows that of the spe- 
cific heat. At high temperatures, where the specific heat 
is constant, A decreases as the phonon mean free path 
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FIG. 3. Thermal conductivity A for a (10, 10) carbon nan- 
otube (solid line), in comparison to a constrained graphite 
monolayer (dash-dotted line), and the basal plane of A A 
graphite (dotted line) at temperatures between 200 K and 
400 K. The inset reproduces the graphite data on an expanded 
scale. The calculated values are compared to the experimen- 
tal data of Refs. |l7j] (open circles), [[L8| (open diamonds) , and 
[p.9[ (open squares) for graphite. 



becomes smaller due to umklapp processes. Our calcu- 
lations suggest that at T = 100 K, carbon nanotubcs 
show an unusually high thermal conductivity value of 
37, 000 W/m-K. This value lies very close to the highest 
value observed in any solid, A = 41, 000 W/m-K, that has 
been reported Q for a 99.9% pure 12 C crystal at 104 K. In 
spite of the decrease of A above 100 K, the room temper- 
ature value of 6,600 W/m-K is still very high, exceeding 
the reported thermal conductivity value of 3, 320 W/m-K 
for nearly isotopically pure diamond 

We found it useful to compare the thermal conductiv- 
ity of a (10, 10) nanotube to that of an isolated graphene 
monolayer as well as bulk graphite. For the graphene 
monolayer, we unrolled the 400- atom large unit cell of 
the (10, 10) nanotube into a plane. The periodically re- 
peated unit cell used in the bulk graphite calculation con- 
tained 720 atoms, arranged in three layers. The results 
of our calculations, presented in Fig. |[ suggest that an 
isolated nanotube shows a very similar thermal trans- 
port behavior as a hypothetical isolated graphene mono- 
layer, in general agreement with available experimental 
data [p7|"^9|. Whereas even larger thermal conductivity 
should be expected for a monolayer than for a nanotube, 
we must consider that unlike the nanotube, a graphene 
monolayer is not self-supporting in vacuum. For all car- 
bon allotropes considered here, we also find that the ther- 
mal conductivity decreases with increasing temperature 
in the range depicted in Fig. [| 

Very interesting is the fact that once graphene lay- 
ers are stacked in graphite, the inter-layer interactions 
quench the thermal conductivity of this system by nearly 
one order of magnitude. For the latter case of crystalline 
graphite, we also found our calculated thermal conductiv- 
ity values to be confirmed by corresponding observations 
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in the basal plane of highest-purity synthetic graphite 
fl7| p^t which are also reproduced in the figure. We would 
like to note that experimental data suggest that the ther- 
mal conductivity in the basal plane of graphite peaks near 
100 K, similar to our nanotube results. 

Based on the above described difference in the con- 
ductivity between a graphene monolayer and graphite, 
we should expect a similar reduction of the thermal 
conductivity when a nanotube is brought into contact 
with other systems. This should occur when nanotubes 
form a bundle or rope, become nested in multi-wall nan- 
otubes, or interact with other nanotubes in the "nan- 
otube mat" of "bucky-paper" and could be verified ex- 
perimentally. Consistent with our conjecture is the low 
value of A«0.7 W/m-K reported for the bulk nanotube 
mat at room temperature [0] . 

In summary, we combined results of equilibrium and 
non-equilibrium molecular dynamics simulations with ac- 
curate carbon potentials to determine the thermal con- 
ductivity A of carbon nanotubes and its dependence on 
temperature. Our results suggest an unusually high value 
A«6,600 W/m-K for an isolated (10,10) nanotube at 
room temperature, comparable to the thermal conduc- 
tivity of a hypothetical isolated graphene monolayer or 
graphite. We believe that these high values of A are as- 
sociated with the large phonon mean free paths in these 
systems. Our numerical data indicate that in presence of 
inter-layer coupling in graphite and related systems, the 
thermal conductivity is reduced significantly to fall into 
the experimentally observed value range. 
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